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Abstract. We report the detection of nonradial g-mode oscillations in the DA white dwarfs HS 1039+4112 (B=15.9) and 
HE 1429-0343 (B=15.8) from time-series photometry made at the Loiano 1.5 m telescope. The two stars were previously 
selected as probable pulsators based on two-color photometry and spectral analysis respectively. Following our temperature and 
surface gravity determinations, HS 1039+41 12 (r cff =l 1200 K, logg=8.2) is located near the red edge of the ZZ Ceti instability 
strip, whereas HE 1429-0343 (T e g=l 1400 K, log g=7.8) falls in the middle of the strip. Both stars show a multi-mode behavior 
with the main periods at about 850 and 970 s respectively, and relatively large amplitudes (~7% and ~2.5%). 
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^ 1 . Introduction of the poorly known rate of the l2 C(a, y) 16 nuclear reaction 

. . ' , „ . , „ „ . . , , s t (see e.g. Metcalfe 2003 and Straniero et al. I2003 and refer- 

►> The ZZ Ceti stars (or DAVs=DA Variables) are old DA (H- 



ences therein), and of the WD cooling rates, including crystal- 
rich) white dwarfs (WDs) showing multi-periodic luminosity ,• . , Rrvrvc j i 

. . inn j 1 onn a nzatlon processes (Kanaan et al. 2005 and ref . therein; see also 



variations with typical periods between ~ 100 and ~ 1200 s and . . , i A A )u i t , . ., , ... ... 

,. , , „ , Fontaine et al. 2001 for a detailed review on this argument). 
r$ amplitudes between a few percent down to less than 0. 1 %, gen- 

' erally interpreted as nonradial g-mode oscillations. ' The effective temperatures of the ZZ Ceti stars are confined 

The importance of studying the ZZ Ceti variables is related between about 10700 and 12100 K for a canonical surface grav- 

to the possibility, given by the seismological analysis, to study '^ }° S S~ 8 (Berger on et al . H223 [2004 ; Koester & Holberg 

in detail their interiors and to derive numerous basic stellar pa- EMD Mukadam et al. 2004a l. The blue edge of the DAV insta- 

rameters such as the mass, the rotation and the thickness of the bilit y stri P can be interpreted as the limit where the base of the 

external layers of H and He. Due to the high stability of the stel- surface convective zone becomes sufficiently deep for the local 

lar structure, the pulsation periods of the DAV white dwarfs are thermal time-scale to be comparable to the shortest observable 

extremely stable in time, with values of P that can be as small «- mode pulsation period (*100 s for typical WD masses, see 

as * 10 15 (Mukadam et al. EOOH and references therein). More Tassoul et al-HSSHand references therein). 

in general, WD seismology can help to improve our knowledge This scenario is in agreement with an observed trend on 

temperature: the hottest DAVs have shorter pulsation periods 

Send offprint requests to: R. Silvotti , , ... , ... , , . „„ 

. , , „ , . and lower amplitudes with respect to the cooler ones. More 

Based on observations collected at the Loiano Observatory, op- .... . . , , , 

... , ,. . „ . ,„,,„ , . .f uncertain is the interpretation of the red edge: convection- 

erated by the Istituto Nazionale di Astrojisica (INAF), and at the r _g_ 

Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, operated Potion interaction ( Winget & Fontaine [TOJ, convective 

jointly by the Max-Planck Institut fur Astronomie and the Institute de mlxln g between hydrogen and helium (Bergeron et al.|1990J or 

Astrofisica de Andalucia (CSIC). Partially based on data obtained at something else ? This is related to the driving of the pulsations 

the Paranal Observatory of the European Southern Observatory for which could be based on the K-y mechanism, active in a hy- 



programmes 165.H-0588 and 167.D-0407 drogen layer near the surface (Dziembowski & Koester u 981 
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Dolez & Vauclair 1981 Winget et al. 119521 . although an al- 
ternative "convective driving" has been proposed by Brickhill 
(119911 see also Gautschy et al. 119961 and Goldreich & Wu 
E999J. 

Recently, thanks to the results of the Sloan Digital Sky 
Survey, the number of known ZZ Ceti pulsators has in- 
creased by a factor of two, reaching a total number of 82 
(Bergeron et al. 120041 Mukadam et al. 2004a Mullally et 
al. |2005 1. Moreover another 5 (possibly 6) DAVs have been 
observed in cataclysmic variable systems (Woudt & Warner 
120041 Warner & Woudt 120051 Araujo-Betancor et al. 120051 . 
Following Mukadam et al. ( 2004b I, it is still not clear whether 
there might be some non-variable stars within the strip, al- 
though these stars could simply have very low pulsation am- 
plitudes, below the detection limit. 

In this article we present two new ZZ Ceti stars, selected 
from the HS (Hamburg Schmidt, Hagen et al. 119951 and HE 
(Hamburg ESO, Christlieb et al. 2001 ) surveys, which bring 
the total number of known non-interacting ZZ Ceti to 84. 

2. Selection of the ZZ Ceti candidates 

The candidates were selected based on two-color photometry 
and spectral analysis. 

2.1. Two-color Stromgren photometry 

A list of 3000 possible cool white dwarf stars was derived from 
the low-resolution photographic objective prism spectra of the 
Hamburg Quasar Survey (HQS) by Homeier & Koester (2001 1. 
One thousand of these objects are cool (T e fft 16000 K) and, be- 
cause of the large errors on the temperature (up to ~6000 K), 
all of them are potential ZZ Ceti pulsators. In order to select the 
best ZZ Ceti candidates, preparatory observations are currently 
being conducted, aimed to derive more precise temperatures 
and gravities for the « 500 brightest stars (B<16.5 mag). 

We obtained photometry in the Stromgren u, b, y and 
Johnson / bands at the Calar Alto 2.2m telescope using the 
BUSCA four-channel CCD camera (Reif et al. IT9991 l . The 
colours u — b and b — y show a strong dependency on T e ff and 
log g, which is strongest in the region of the ZZ Ceti instabil- 
ity strip. In the course of our analysis, we noticed technical 
problems with the BUSCA y filter as well as an even stronger 
temperature dependency of b - I than that of b - y. Therefore 
we use the u — b vs. b — I two-color diagram for the analysis of 
our data. 

To derive atmospheric parameters from the two-color data, 
we compare them to a grid of synthetic colors. This grid was 
computed from Koester DA model atmospheres (described e.g. 
in Finley et al. 1997) by convolving the model spectra with 
the BUSCA filter transmission curves. Temperatures and gravi- 
ties were determined by interpolating linearly between the grid 
points. The observational and calibration errors of the colors 
are ~ 0.02 mag, which translates into an error of ~ +0.1 in 
logg and ~ +300 K in T e ff. This temperature error is small 

1 See also the web pages of BUSCA at 
http://www. caha. es/guijarro/BUSCA/intro. html 
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Fig. 1. Two-color data and the grid of synthetic colors from 
model atmospheres. Solid lines are lines of constant gravity and 
dotted lines are lines of constant temperature (r e jf=20000 K, 
18000, 16000, 15000 and down to 9000 in steps of 1000 K). 
Dashed lines show the boundaries of the instability strip, 
according to Bergeron et al. (2004). The crosses represent 
the four stars selected as possible pulsators: from left to 
right HS 1610+1639, HS 1443+2934, HS 1253+1033 and 
HS 1039+41 12. The squares are two previously known ZZ Ceti 
stars, GD 244 (left) and PG 2303+242, that were included in 
the preparatory observations as a check of our methods. 

compared to the width of the instability strip, and thus most 
objects are clearly placed either outside or inside the instability 
strip by the photometry. 

Four objects that were selected as potential pulsators in 
this way are studied here. Their two-color data are plotted in 
Fig. 1, together with the temperature/gravity grid and the in- 
stability strip edges, taken from Bergeron et al. (2004). For 
one of the new pulsators, HS 1039+41 12 (RA 20 oo=10 42 33.5, 
DEC 2 ooo=+40 57 16, B= 15. 9 mag), we find r eff =l 1200+270 K 
and log g=8.2±0. 12. The parameters of the three remaining ob- 
jects are given in Table 2. 

The selection of the object HS 1610+1639 is problematic: 
using the now obsolete b - y data, we obtained r e ff=l 1700 K 
and log g=8.2, corresponding to the middle of the instability 
strip. This is the reason why HS 1610+1639 was included in 
the sample observed at Loiano. However, after the time series 
observations at Loiano had been made, the re-analysis of the 
data, now using the b-I color, has yielded a higher temperature 
for this object, T eB - 14500 K (and logg=8.4), definitively too 
hot for ZZ Ceti pulsations. 

2.2. Spectral analysis 

The other new pulsator, HE 1429-0343 (RA 20 oo=14 32 03.2, 
DEC2ooo=-03 56 38, B=15.8), was selected based on high- 
resolution spectra from the SPY survey (Napiwotzki et al. 
2003). An analysis of these spectra for 754 white dwarfs, many 
of which are first spectroscopic confirmations of WDs from the 
HS and HE surveys, was conducted by Voss et al. (2005, in 
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preparation). Details about theoretical spectra and fitting pro- 
cedures can be found in Finley et al. (1997). Fig. 2 displays 
the model atmosphere fit for HE 1429-0343, from which we 
obtain r eff =l 1434 + 36 K and log g=7.82 + 0.02. 
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Fig. 2. The SPY spectrum and fit of HE 1429-0343. The spec- 
trum was rebinned to a resolution of 1A for this plot, the fit 
was done with the original resolution of 0.1 A. The deep core 
of Hy6 is not real but is caused by a cosmic, which however does 
not significantly affect the fitting process. 
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Fig. 3. Spectrum of HS 1641+7132 and fits for T eff = 15 350 K, 
log g = 7.9 (solid line), and for 1 1 800 K, log g = 8. 1 (dashed) 



The other object selected from spectroscopy is 
HS 1641+7132, one of the 12 possible DAVs identified 
from medium-resolution spectroscopy by Homeier et al. 
( 1998 1. HS 1641+7132 was selected for time-series photome- 
try because of its atmospheric parameters, T e g=l 1600 ± 80 K 
and logg=8.0 ± 0.06, which correspond to the middle of the 
instability strip. While the formal errors are small despite the 
low S IN and the resolution (6 A), line profile fits near the 
Balmer line maximum are known to easily lead to degenerate 
solutions for T e s and logg. Indeed, upon re-analyzing the 



original data from 119981 we found that even small changes 
to the startup values of the fitting procedure would produce 
temperatures as much as 200 K apart, with a mean value of 
1 1800 K. In addition, more important, we found a second solu- 
tion for the fit at much higher temperature T^- 15350 ± 300 K 
and log g=7.88 ± 0.06, with a reduced ^2= 1.26 as compared 
to 1 .45 for the cooler solution. The two fits of the spectrum are 
shown in Fig. 3. 

3. Time-series photometry 

The time-series observations were performed in June 2004 at 
the 1.5 m Loiano telescope using the BFOSC 2 instrument 
without any filter. A recent version of the software allows to 
launch a series of exposures where the CCD is read only in a 
small box around the target, in order to decrease the read-out 
time. With a box of typically ~300x400 pixels, we were able to 
measure at least two reference stars in each field. The integra- 
tion times used were between 12 and 18 s, giving an effective 
resolution of about 24 to 30 s, because of the CCD reading plus 
some extra time due to internal controls of the software. 

The data were then reduced using IRAF 3 aperture pho- 
tometry. The flux of each target was divided by a combination 
of the reference stars with different weights, and then normal- 
ized to the mean ratio. Finally a residual extinction correction 
was applied and the times were converted to Barycentric Julian 
Date (BJD), using the algorithm of Stumpff ( 1980). 

The light curves of the two new ZZ Ceti stars are shown 
in Fig. 4 and 5. The variable amplitudes in the light curves, 
presumably due to beating effects, suggest that both targets 
have more than one excited pulsation mode. Each target was 
observed in two different nights to confirm the pulsations. In 
this way, joining together the two runs, the frequency resolu- 
tion of the Fourier Transform is improved and it is possible to 
make a first separation of the different peaks (Fig. 4 and 5). The 
light curve of HS 1039+4112 has the typical large amplitude 
of many of the cool ZZ Ceti; moreover the spectrum shows the 
possible presence of the first and second harmonic of the main 
peak at 1.181 mHz (846.6 s). In the case of HE 1429-0343, 
which is located more close to the center of the instability strip, 
the pulsation amplitude is lower. Close to the highest peak of 
the spectrum at 1.032 mHz (969.0 s) 4 there are at least two 
other signals at 0.922 and 1.206 mHz (1084.9 and 829.3 s re- 
spectively). 

In Table 1 we report an attempt to fit the observed light 
curves (both nights together) with sinusoidal fits; the results 
of the fits are also shown in the upper panels of Fig. 4 and 5. 
However, due to the poor coverage, the pulsation spectra are 
not resolved and therefore these results should be taken with 
some prudence, in particular regarding the frequencies with 

2 Bologna Faint Object Spectrograph & Camera, see the web site 
of the Loiano Observatory at http://www.bo.astro.it/loiano/ for more 
details. 

3 Image Reduction and Analysis Facility, written and supported 
by the IRAF programming group at the National Optical Astronomy 
Observatories (NOAO) in Tucson, Arizona l http://iraf.noao.edu/ 1. 

4 The amplitude of this mode, however, is not the highest following 
the results of the sinusoidal fits, as reported in Table 1. 
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Fig. 4. Light curves, spectral window and amplitude spectrum 
of HS 1039+4112 (1 mma = 0.1% change in intensity). The 
observations were performed in June 20 and 22, 2004. The am- 
plitude spectrum is obtained joining together the two runs in 
order to increase the frequency resolution. The lower panels 
show the behavior of the spectrum when we subtract, at each 
iteration, the highest peak (pre-whitening); note that the inten- 
sity scale is not the same in all panels. 



lower amplitudes (3F1 and F4), which are quite close to the 
noise level. Concerning HS 1039+4112, the first and second 
harmonics suggest that the true value of the main frequency 
could be 1.169 mHz, whereas the value of 1.181 detected in 
the power spectrum could correspond to the 1-day alias (1.169 
+ 0.012 mHz). Indeed, using 1.169 mHz, the least-squares fit is 
more stable and therefore we adopt this solution in Table 1 (the 
alternative solution is also given in square brackets for com- 
pleteness). In the case of HE 1429-0343, the results of the si- 
nusoidal fits favor a higher amplitude of F3 and Fl respect to 
F2 which, however, is the highest peak in the Fourier spectrum 
(24.9 mma). 



Fig. 5. Same as Fig. 4, but for HE 1429-0343. The observa- 
tions of HE 1429-0343 were performed in June 18 and 21, 
2004. 



Table 1. Results of the non-linear least-squares fits 



Name 




F(mHz) 


P(s) 


A(mma) 


HS 1039+4112 


Fl 


1.169 


855.5 


55.2 






[1.181] 


[846.6] 


[71.8] 




F2 


1.194 


837.3 


26.0 










[30.7] 




2F1 


2.338 


427.8 


29.1 




3F1 


3.507 


285.1 


14.6 


HE 1429-0343 


Fl 


0.922 


1084.9 


16.3 




F2 


1.032 


969.0 


12.7 




F3 


1.206 


829.3 


18.3 




F4 


2.222 


450.1 


10.2 
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4. The ZZ Ceti instability strip updated 

In Fig. 6 the location of the two new ZZ Ceti stars in a (7^, 
log g) plane is shown and compared with the edges of the obser- 
vational ZZ Ceti instability strip as defined by Bergeron et al. 
( 2004 1. Both stars are within the strip, although HS 1039+41 12 
is very close to the red edge. 
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Fig. 6. The observational ZZ Ceti instability strip updated: the 
two new pulsators are represented by the two asterisks with 
error bars, whereas the two objects close to the blue edge 
are HS 1253 + 1033 and HS 1443+2934 (from left to right), 
for which no pulsations were found. The filled circles are 
the known ZZ Ceti stars from Bergeron et al. (2004), the 
open circles are the new pulsators discovered from the SDSS 
(Mukadam et al. 2004a; Mullally et al. 2005). 

On the other hand, if we consider the other four targets ob- 
served in the same run at Loiano, none of them shows pulsa- 
tions. Their upper limits (i.e. the amplitude of the highest peak 
in the region of interest of their temporal spectra) are reported 
in Table 2, together with their B magnitudes, effective temper- 
atures and surface gravities, date and duration of the observa- 
tions. 

Concerning HS 1641+7132, its non-variability is compati- 
ble with the higher temperature solution which was found upon 
re-analyzing the spectroscopic data, as described in section 
2.2. Therefore we conclude that this star is likely to be a non- 
variable DA with T eB between 15000 and 16000K. The same 
applies to HS 1610+1639, which is clearly far from the blue 
edge of the instability strip, having an effective temperature of 
about 14500 K. 

Finally, for what concerns the other two stars which 
do not show luminosity variations, HS 1253+1033 and 
HS 1443+2934, both lie near the blue edge. The case of 
HS 1253+1033 is the most intriguing as this star apparently 
falls inside the observational strip of Bergeron et al. (|2004i. 



However, even in this case, the relatively large uncertainties 
in r e ff and logg can not completely exclude the possibil- 
ity that HS 1253+1033 is actually slightly beyond the blue 
edge. Moreover, looking at Fig. 6, both stars would be outside 
the observational strip as defined by the SDSS sample alone. 
Considering that we used the same Koester model atmospheres 
as used by the SDSS team, the stability against g-modes of 
HS 1253+1033 and HS 1443+2934 appears still less surpris- 
ing. 
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Table 2. Objects which do not show pulsations 



Name 


RA2000 


DEC2000 


B 




l°g£ 


UT date 


D 


L 










(K) 


(cgs) 


at start 


(hours) 


(mma) 


HS 1253+1033 


12 56 28.5 


+10 17 09 


14.4 


12600 


8.5 


16/06/04 


1.6 


5.3 


HS 1443+2934 


14 45 28.5 


+29 21 31 


14.5 


12400 


8.1 


22/06/04 


2.0 


1.5 














23/06/04 


1.5 


1.2 


HS 1610+1639 


16 13 02.4 


+ 1631 56 


15.8 


14500 


8.4 


17/06/04 


3.4 


2.0 


HS 1641+7132 


16 41 00.0 


+71 26 59 


16.7 


15300* 


7.9* 


23/06/04 


1.4 


4.6 



*An alternative solution for this star gives r efr = 11600 K and logg=8.0 (see section 2.2 for more details). 



Winget D.E., Fontaine G., 1982, in Pulsations in Classical and 
Cataclysmic Variable Stars, ed. J.P. Cox and J. Hansen (Univ. of 
Colorado: Boulder), p.46 

Winget D.E., Van Horn H.M., Tassoul M„ et al., 1982, ApJ 252, L65 

Woudt P.A., Warner B., 2004, MNRAS 348, 599 



